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Abstract. We suggest a novel picture of the quantum Universe — its creation is 
described by the density matrix defined by the Euclidean path integral. This yields 
an ensemble of universes — a cosmological landscape — in a mixed state which 
is shown to be dynamically more preferable than the pure quantum state of the 
Hartle-Hawking type. The latter is dynamically suppressed by the infinitely large 
positive action of its instanton, generated by the conformal anomaly of quantum fields 
within the cosmological bootstrap (the self-consistent back reaction of hot matter). 
This bootstrap suggests a solution to the problem of boundedness of the on-shell 
cosmological action and eliminates the infrared catastrophe of small cosmological 
constant in Euclidean quantum gravity. The cosmological landscape turns out to be 
limited to a bounded range of the cosmological constant A m j n < A < A max . The domain 
A < A m ; n is ruled out by the back reaction effect which we analyze by solving effective 
Euclidean equations of motion. The upper cutoff is enforced by the quantum effects 
of vacuum energy and the conformal anomaly mediated by a special ghost-avoidance 
renormalization of the effective action. They establish a new quantum scale A max which 
is determined by the coefficient of the topological Gauss-Bonnet term in the conformal 
anomaly. This scale is realized as the upper bound — the limiting point of an infinite 
sequence of garland-type instantons which constitute the full cosmological landscape. 
The dependence of the cosmological constant range on particle phenomenology suggests 
a possible dynamical selection mechanism for the landscape of string vacua. 
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1. Introduction 

The ideas of quantum cosmology [H [2] and Euclidean quantum gravity [SI Ej are again 
attracting attention. One of the reasons is the fact that the landscape of string vacua is 
too big jH] to hope that a reasonable selection mechanism for the plethora of these vacua 
can be successfully worked out within string theory itself. Thus, it is expected that other 
methods have to be invoked, at least some of them appealing to the construction of the 
cosmological wavefunction [HI El E] • 

This approach is based on the idea of quantum tunneling from the classically 
forbidden state of the gravitational field. Semiclassically this state is described in 
terms of the imaginary time, that is by means of the Euclidean spacetime, so that the 
corresponding amplitudes and probabilities are weighted by the exponentiated Euclidean 
gravitational action, exp(— S*e). The action is calculated on the gravitational instanton 
- the saddle point of an underlying path integral over Euclidean 4-geometries. This 
instanton gives rise to the Lorentzian signature spacetime by analytic continuation 
across minimal hypersurfaces, and the instanton topology allows one to interpret this 
continuation either as a quantum tunneling or creation of the Universe from "nothing" . 
Then the minima of Se give most probable values of the physically interesting parameters 
of the Universe, including initial conditions and fundamental coupling constants [TUl lllj. 
Apparently, when applied to the string landscape this framework can serve as a selection 
mechanism for the most probable physical vacua. 

There is an immediate difficulty with this program, which is rooted in the problem of 
unboundedness of the Einstein action. The most important example of the creation from 
"nothing" is the Hartle-Hawking wave function of the Universe 0111, which describes 
nucleation of the de Sitter Universe from the Euclidean 4-dimensional hemisphere, 



It has a negative action which diverges to — oo for the cosmological constant A — > 0. 
This means that the vanishing cosmological constant is infinitely more probable than a 
positive one - the result which has a very controversial status. On one hand it was used 
as a justification of a zero value of A [10J by the Coleman big-fix mechanism [H], but 
on the other hand it was considered as extremely confusing and anti-intuitive because 
it preferred creation of infinitely large universes and disfavored inflation. No need to 
say that this result also looks unappropriate in context of the cosmological acceleration 
which requires very small, but still nonvanishing cosmological constant (certainly, if we 
identify dark energy phenomenon with the cosmological constant effect). 

Such a controversy gave rise to the so-called cosmology debate [12] reflecting the 
first attempt to circumvent this problem. Point is that the alternative formalism of 
the Wheeler-DeWitt equation, which is only formally equivalent to the path integral 
approach [TSJ H3|, admits also an opposite sign in the exponential of (JTJ). It exists 
either in the form of Linde ^H] or Vilenkin ^B] tunneling proposals. These proposals 
differ by boundary conditions ^2] but qualitatively render the same exponentially small 
probability for A — > 0, which opens the possibility for conclusions opposite to the Hartle- 




(1) 
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Hawking case. In particular, the inclusion of one-loop effects allows one to shift most 
probable values of the effective cosmological constant from zero to a narrow highly 
peaked range compatible with the requirements of inflation [TT| . 

Such a big discrepancy in handling this infrared domain gives a feeling of 
indeterminacy - fundamentals of quantum gravitational tunneling seem to be on too 
shaky a ground to serve as a reliable completion of the string theory of everything. 
Apparently, this discrepancy could be resolved within a deeper understanding of the 
gravitational path integral and correct operator realization of the Wheeler-DeWitt 
equation. However, this achievement might still be in vain, because in view of the 
anticipated string nature of everything this equation might be not fundamental. For 
this and other reasons (we will not repeat here usual arguments refuting the tunneling 
prescription on the basis of erroneousness of anti-Wick rotation, etc.) we will stick to 
the Hartle-Hawking prescription of the Euclidean path integration. 

In this prescription the attempts to resolve the infrared catastrophe of small A 
are very hard to implement. Indeed, within a conventional wisdom of a low-derivative 
expansion all quantum corrections are of higher order in the curvature R than the 
cosmological and Einstein terms. Therefore, they are negligibly small in the infrared 
limit R ~ A — * and cannot compete with the infinitely growing tree-level contribution 
of This is a very general argument that cannot be circumvented within the local 
curvature expansion of the effective action. So the only reasonable option is to try 
essentially nonlocal quantum effects mediated by nonlocal terms of non-vacuum nature. 

This idea was recently attempted in speculative but very thought-provoking papers 
of Henry Tye et al [7J E] who suggested to consider the effect of radiation in quantum 
creation from nothing (see also [9J). Radiation effectively appeared in [?J E] as a result of 
ultraviolet (and infrared) renormalization of the effective action of the "environment" - 
spatially inhomogeneous modes of matter and gravitational fields, treated perturbatively 
on the minisuperspace FRW background. Their contribution to the instanton action was 
advocated in [7JIE] to be positive and inverse proportional to the square of A, so that the 
action was bounded from below at some positive A with which the Universe gets created 
from "nothing". This mechanism has already been justifiably criticized in [TH] by noting 
that this type of "radiation" can only renormalize the effective value of A, in terms of 
which the boundedness properties of the action remain the same. Nevertheless, as we 
show below, the radiation contents of the creation-from-nothing instanton can really 
mediate the boundedness of its effective action. However, the origin of this radiation 
and the mechanism of the infrared cutoff at small A are drastically different from those 

of [HE]. 

The core of our suggestion is a simple observation that the presence of radiation 
implies a statistical ensemble described by the density matrix, rather than a pure state 
assumed in [TIE]- Density matrix in Euclidean quantum gravity [THj, p[p,p'], originates 
from an instanton with two disjoint boundaries E and E' associated respectively with 
its two entries ip and <p' (collecting both gravitational and matter variables). The 
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Figure 1. Picture of instanton representing the density matrix. Dashed lines 
depict the Lorentzian Universe nucleating from the instanton at the minimal 
surfaces £ and 



bulk part of the instanton interpolates between these boundaries and, thus, establishes 
mixing correlations between the observables located at them, see Fig^ In contrast, the 




Figure 2. Density matrix of a pure state represented by the union of two 
vacuum Hartle-Hawking instantons giving rise to Lorentzian spacetime. 

pure density matrix of the Hartle-Hawking state corresponds to the situation when the 
instanton bridge between E and E' is broken, so that topologically the instanton is a 
union of two disjoint hemispheres. Each of the half-instantons smoothly closes up at 
its pole which is a regular internal point of the Euclidean spacetime ball (Fig|2I) — a 
picture illustrating the factorization of p = I^hhX^hhI- 

The concept of an initial state of the Universe should be extended to include mixed 
states of the above type. This extension seems natural because it is up to Nature 
to decide whether a pure or mixed initial state is preferable. So this is a dynamical 
question rather than a postulate, and we address this question below by incorporating 
both types of states into a unified framework of the density matrix, t. This is attained by 
considering the density matrix as a Euclidean path integral over the fields on spacetimes 

X It should be emphasized that, unlike in black hole theory and the deSitter cosmology with horizons, 
mixed nature of the quantum state is fundamental here, because it is not effectively acquired as a result 
of tracing out invisible superhorizon degrees of freedom or coarse graining. This is just a fundamental 
lack of physical information about the quantum state which should be described in terms of a density 
matrix rather than a pure quantum state. This situation is not postulated but treated on equal footing 
with the case of a pure quantum state carrying the full amount of physics information admissible by the 
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with a tubular topology of FigJTJ which includes as a limiting case the topology of Fig|2] 
with a broken bridge. 

Qualitatively, the main effect that we advocate here is as follows. When 
calculated in the saddle-point approximation the density matrix automatically gives 
rise to radiation whose thermal fluctuations destroy the Hartle-Hawking instanton. 
Namely, the radiation stress tensor prevents the half-instantons of the above type 
from closing and, thus, forces the tubular structure on the full instanton supporting 
the thermodynamical nature of the physical state. The existence of radiation, in 
its turn, naturally follows from the partition function of this state. The partition 
function originates from integrating out the field <p in the coincidence limit tp' = p. 
This corresponds to the identification of E' and E, so that the underlying instanton 
acquires toroidal topology, see FigHU Its points are labeled by the periodically 
identified Euclidean time, a period being related to the inverse temperature of the 
quasi-equilibrium radiation. The back reaction of this radiation supports the instanton 
geometry in which this radiation exists, and we derive the equation which makes this 
bootstrap consistent. 




E 

Figure 3. Calculation of the partition function represented by compactification 
of the instanton to a torus with periodically identified Euclidean time. 

We show that for the radiation of conformally-invariant fields the analytical and 
numerical solution of the bootstrap equations yields a set of instantons — a landscape 
— only in the bounded range of A, 

A min < A < A max . (2) 

This set consists of the countable sequence of one-parameter families of instantons 
which we call garlands, labeled by the number k = 1,2,3,... of their elementary 
links (of the form shown on FigGJ). Each of these families spans a continuous subset, 
A^ n < A < Ama X , belonging to (J2J). These subsets of monotonically decreasing lengths 
Ama, x — A^ m ~ l/k 4 do not overlap and their sequence has an accumulation point at 
the upper boundary A ma x of the range (0). Each of the instanton families at its upper 

principles of quantum mechanics. As we will see below, the system dynamically prefers mixed states 
of quasi-thermal nature. 
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boundary Amax saturates with the static Einstein Universe filled by a hot equilibrium 
radiation with the temperature T(&) ~ mp/ Ink 2 , k 3> 1, and having the negative 
decreasing with k action 1^ ~ — \n 3 k 2 /k 2 . This means that large k garlands do not 
dominate the cosmological landscape and turn out to be less probable than k = 1 
instantons. 

Remarkably, all values of A below the range (J2J), A < A min , are completely ruled 
out either because of the absence of instanton solutions or because of their infinitely 
large positive action. This effect is provided by the radiation back reaction and survives 
even in the classical limit of the theory. A similar situation holds for A > A max - no 
instantons exist there, and the Lorentzian configurations in this overbarrier domain (if 
any) are exponentially suppressed relative to those of (j2J. However, the high-energy 
cutoff is entirely the quantum effect of the vacuum energy of quantized fields and their 
conformal anomaly. This effect generates a new scale A max in gravity theory. Specifically, 
this scale is determined by the coefficient of the (topological) Gauss-Bonnet term in the 
conformal anomaly. 

Thus, neither infrared nor ultraviolet catastrophe is present in our picture. Rather, 
we get a bounded cosmological landscape. Its boundaries depend on the phenomenology 
of quantum fields. In particular, for a single conformal scalar and vector field they are 
respectively given by the following Planckian values, nip = 37r/4G, 

A£f« 8.991 m%, A^ ar = 360m|, (3) 
A^ n tor « 3.189 m 2 P , A^r = 180/31 m 2 P . (4) 

Qualitatively, these bounds decrease with the growth of the field spin and also go to 
zero as 1/N with the growing number of fields N. This brings us to the domain of 
semilassical consistency of the model. This also means that the landscape shrinks to a 
narrower and narrower subplanckian range when we ascend the phenomenology energy 
scale. This suggests the possibility of generating a long sought selection mechanism 
which might disentangle from the vast string landscape a dynamically preferred vacuum 
compatible with the observed Universe. 

The paper is organized as follows. In Sect. 2 we present the effect of radiation in 
Euclidean cosmology which maintains a thermal nature of the cosmological instanton 
and underlies the density-matrix entanglement. In Sect. 3 we construct the path integral 
representation for this density matrix and the corresponding partition function (or 
the Euclidean effective action). Here we also describe the approximation when the 
minisuperspace metric of the instanton is subject to the self-consistent back reaction of 
quantum matter. In Sects. 4-6 we calculate the effective action of conformally invariant 
fields on the instanton with generic FRW metric and derive the nonlocal effective 
equation for the latter. The effective action is obtained as a sum of contributions of the 
conformal anomaly of quantum fields and their action on the static Einstein Universe. 
They are both properly renormalized to avoid ghost instabilities in the solutions of 
the corresponding effective equation. The latter turns out to be the Friedmann 
equation which is modified by the conformal anomaly term and has as a source the 
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radiation density term determined by the bootstrap equation. Elimination of the 
infrared catastrophe by the back-reaction effect of a conformal matter is demonstrated 
within this bootstrap in Sect. 7. The infinite sequence of garlands - multiply folded 
instantons generating the new quantum scale as an upper bound of the landscape (J2J) - 
is constructed in Sect. 8. The concluding section summarizes the results and discusses 
the extension of the model to quantum fields of higher spins and a possible application 
of the bounded cosmological landscape as a means of reducing the landscape of string 
vacua. 



2. The effect of radiation 



The effect of radiation on cosmological solutions of Euclidean gravity is as follows. In a 
spatially closed cosmology with the Euclidean FRW metric (f^ 3 -* is a 3-sphere of a unit 
radius) 

als 2 = N 2 (r)dr 2 + a 2 (r)d 2 n^, (5) 

the Friedmann equation includes the radiation density term inverse proportional to the 
fourth power of the scale factor a 

- 2 = ^~H^-,. (6) 

or cr ar 

Here we choose the gauge of a unit lapse N — 1, H 2 = A/3 is the cosmological constant 
in terms of the Hubble parameter H and the constant C characterizes the amount 
of radiation. The radiation term blows up to infinity for a — > and in view of the 
positivity of the radiation density (C > 0) prevents the model from the cosmological 
collapse. Moreover, the right hand side here also vanishes for some big a, and the 
Friedmann equation has two turning points - the situation analyzed in the context of 
quantum tunneling in ffiH I2T] . In the corresponding solution 



1 



a(r) = -^\J l ~ C 1 -4C# 2 )V2 CO s 2Hr (7) 
the scale factor never goes to zero and varies between the maximal and minimal values 



a ± = -jLjy/l±(l-4CIP)V* (8) 

during one period of the Euclidean time —it/2H < t < tt/2H. This instanton solution 
exists for a finite amount of radiation satisfying the bound 

AH 2 C < 1, (9) 

and it can be used for constructing the density matrix of the Universe if we identify the 
boundaries £ and £' with the minimal surfaces of the maximal size a + at r = ±tt/2H. 
The Lorentzian spacetime with the quasi-deSitter metric nucleates from them via 
the analytic continuation in time and gives rise to inflationary cosmology with the 
asymptotic value of the Hubble constant H at late stages of expansion, see FigJT] above. 
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As we see, at any nonzero H the radiation term prevents the instanton bridge between E 
and £' from breaking and thus supports the thermal nature of the ensemble of universes. 

Note that the (half-period) instanton interpolating between the minimal and 
maximal values of a scale factor, which was used in [201 E] for the description of 
quantum tunneling between different Lorentzian domains, is not suitable for the density 
matrix. The calculation of thermodynamic averages incorporates taking the trace and 
leads to the identification of £ and £' surfaces. Therefore, these surfaces should have the 
same geometric characteristics. Thus, only full-period and multiple-period instantons 
are admissible. 



3. The density matrix 

The picture of the above type serves as a tree-level approximation for the density matrix 
of this ensemble. At the exact level this density matrix is given by the Euclidean path 
integral 

P W^']=e T J L>[s,0] exp (-5^,0]). (10) 

9,<t>\ £,£' = ( VW') 

Here Se[<7, <P\ is the classical action of the theory, and the integration runs over 
gravitational g and matter <ft fields on a tubular spacetime interpolating between the 
field configurations ip and ip' at £ and £' (which, as mentioned above, also include both 
gravitational and matter variables). The normalization factor exp(r) is determined from 
the condition trp = 1. This condition serves as a definition of the Euclidean effective 
action 

e~ r = J L>[<7,0]exp(-£ E [<7,0]) (11) 

g,4>\ s = g,<t>\x> 

on the compactified spacetime obtained by identifying the boundaries £ and £' as shown 
on Fig 01 above. Now the integration runs over periodically identified gravity and matter 
fields. 

The motivation for this definition of the density matrix and its statistical sum can 
be found in ^H]- Here we would only mention that this is a natural generalization 
of the path integral for the no-boundary wave function of the Universe. Also, it can 
be regarded as a gravitational generalization of the density matrix of the equilibrium 
thermodynamical ensemble at finite temperature T = 1//3, p = exp (r — (3H), for 
a system with a Hamiltonian operator H. Its kernel in the functional coordinate 
representation p[<p, <p'] = (<p\ p\<p'), similarly to (fTU|) . is given by the Euclidean path 
integral over histories 0(r) in the imaginary time r, interpolating between 0(0) = tp' 
and 4>(/3) = (p. 

This definition suggests an unusual interpretation of the cosmological ensemble. 
It follows that, semiclassically, the set of Lorentzian spacetime universes originates by 
analytic continuation into the complex plane of the temperature variable rather than 
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the time variable. Therefore, the whole picture looks more like a thermal state of an 
already existing universe, rather than its quantum birth from "nothing" by means of 
tunneling from the classically forbidden state. In fact, this situation is a manifestation 
of the principle that "the mathematical formalism [of the quantum theory] is capable of 
yielding its own interpretation" put forward by DeWitt in the other context [22J. Indeed, 
one and the same formalism of the imaginary time is used to describe qualitatively 
different phenomena - - thermodynamics and quantum tunneling. On the physical 
grounds they are easily distinguishable in usual field-theoretic problems, but here they 
exhibit intrinsic dualism between thermodynamics and quantum tunneling. § 

Anyway, the density matrix (fTT)|) unambiguously prescribes a particular mixed 
quantum state of the system. Its semiclassical calculation yields as a saddle point the 
configuration of a tubular topology depicted on FigHJ Its Euclidean part (interpreted 
either in thermodynamical or tunneling terms) is bounded by two minimal surfaces 
S and £'. The analytic continuation across these surfaces yields the ensemble of 
cosmological models expanding in real Lorentzian time, and this picture can be called 
as the origin of cosmological thermodynamics via creation from "nothing". Taking 
the trace of p in the normalization condition for the density matrix results in the 
identification of £ and £' and the toroidal compactification of the instanton depicted 
on FigEl This underlies the semiclassical calculation of the statistical sum and the 
corresponding Euclidean effective action pijl . 

The thermal quasi-equilibrium nature of this statistical sum will be demonstrated 
by explicit calculations in Sect. 5 below. However, the thermodynamical analogy for 
the density matrix kernel, that was mentioned above, is not entirely precise here due 
to general covariance of gravity theory. This happens because: i) the Hamiltonian of a 
spatially closed cosmology is vanishing and ii) the temperature in (fTT)|) is not specified by 
hands as a free parameter, but rather integrated over (as an integration over the lapse 
of Euclidean time separating the surfaces £ and £'). Within semiclassical expansion the 
second property works as follows. The effective "temperature" of the instanton (or its 
Euclidean time period) is uniquely determined from its equation of motion - Euclidean 
Friedmann equation for a cosmological scale factor. 

This equation and the way it incorporates back reaction of quantum matter arise 
as follows. We decompose the full configuration space of [g,4>] in (JUJ), [g,<f>] — > 
[g (T),§(x)}, into a minisuperspace of the cosmological scale factor a(r) and the 
lapse function N(r), go(r) = (a(r), N(r)), and the sector of "matter" $(x) which 
includes together with matter fields also the metric perturbations on the minisuperspace 
background (jSJ), $(x) = (<fi(x), ip(x), A^x), h^ u (x), ...). Under this decomposition the 
measure factorizes accordingly D[g,(f>] = Dg^r) x D§(x) (both factors implying 
relevant gauge-fixing procedures in minisuperspace and field sectors respectively), and 

§ Similarly, the formalism of quantum black holes can be interpreted in terms of widely accepted black 
hole thermodynamics and much less recognized physics of virtual black holes 23 . 
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the statistical sum (|11|) takes the form 

e~ r = J D 9o (t) exp (-r[<?o(T)]), (12) 

where the effective action r[(7o( r )] of quantized matter on the minisuperspace 
background go{r) is determined by the path integral 

e -r[g (r)} = J m{x) exp(-5 E [( ?0 (r),$(x)]). (13) 

Our approximation in what follows will be to include into this action the one-loop 
order, 

r[^] = 5 , E [</] + ri_i oa p[(7], (14) 

and calculate the minisuperspace integral (fT2*|) in the tree-level approximation, which 
reduces to solving the effective equations for r [ go (r) ] . This will give the lowest order 
back reaction effect of quantum matter in the process of quantum creation of the 
cosmological ensemble. 



4. Conformal anomaly and ghosts 

Even for a simple FRW background the one-loop action is not exactly calculable for 
quantum fields of a general type. For conformally-invariant fields, however, one can 
apply the technique of the conformal transformation which relates the FRW metric (J5|) 
rewritten in terms of the conformal time n, 

als 2 = a 2 (7])(dr] 2 + d 2 n^), (15) 

to the metric of the Einstein static Universe of a unit size 

ds 2 = dri 2 + d 2 n®. (16) 

In contrast to the local conformal invariance of the classical action of a conformal 
field of weight w, S conf [e u g^(x), e~ WCT / 2 <p{x) ] = S conf [g^(x), <f>[x) ], a = a{x), its 
quantum effective action has a conformal anomaly 

gJ 1 ^ = ^9 1/2 (aDR + ( 3E + jC 2 ^) , (17) 

determined by the coefficients of the scalar OR, Weyl tensor squared term C 2 ua/3 and 
the Gauss-Bonnet invariant 

E = R lua 7 ~ ^ R lu + r2 - ( 18 ) 

These coefficients depend on the spin of a conformal field j2Hj and read respectively for 
a scalar, Weyl spinor and vector field as 

18 124 
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where for the vector case we cite the result of the ^-function [2^ and point separation 
|2*7] regularizations, which differs from the dimensional regularization result [2H] for a. 

This anomaly, when integrated functionally along the orbit of the conformal group, 
gives the relation between the actions on generic conformally related backgrounds j25] 

r^ loop [g} + AT[g,g}, g^{x) = e a ^g^{x), (20) 



where |3Ti] 



Fi-ioo P [g I 



^[9,9] 



2(4tt) 



dW /2 \ ~ 



2 

liua/3 



+ P[E- 2 -UR 



a 



2 



{Uaf + 2 -R{V,af 



a 
12 



+ 



P 



fd'x (g^R\g)-g^R\g)), {21) 



2(4tt) 2 V12 ' U 

and all barred quantities are calculated with respect to the metric g^ u (x). 

For the conformal factor e CT = a 2 (r) this expression immediately generates higher- 
order derivative terms ~ a 2 in the effective Lagrangian. This certainly produces ghost 
instabilities in solutions of effective equations. It is remarkable, however, that higher- 
derivative terms are all proportional to the coefficient a, because similar terms linear 
in /3 completely cancel out in ()21|) ((Dcr) 2 against the _R 2 (g)-term). The a-term of the 
conformal anomaly can be arbitrarily changed by adding a local counterterm 
admissible from the viewpoint of general renormalization theory. Therefore we can put 



it identically to zero, 
the effective action 



a 



cir = 0, by the following ghost- avoidance renormalization of 



ri-ioop[^]^r^[^] = r 1 _ loop [( 7 ] + ^-p^ J d'xg 1 ' 2 R 2 ( g ). (22) 

Certainly, this additionally spoils conformal invariance of the theory which was anyway 
irreversibly broken by quantum corrections. Thus it is reasonable to fix this local 
renormalization ambiguity by the additional criterion of the absence of ghosts, what 
we do here for sake of consistency of the theory at the quantum level. || 

From (|21|) it then follows that the conformal contribution to the renormalized action 
on the minisuperspace background (fTB]) - (fTBj) equals 

^ 2 1 a 4 ' 



AT[g,g} = T R [g 
3 



r 



Rl9 



m. 



B / dr 



a 
a 



6 a 



mfpB 



(23) 
(24) 



with the constant m 2 p B which for scalars, two-component spinors and vectors equals 
respectively 1/240, 11/480 and 31/120. As we will see, the positivity of the first term 

|| This is certainly not an exhaustive solution of the ghost problem in effective equations of local 
field models. Higher-derivative terms still remain in the other sectors of the theory - the graviton 
sector of transverse-traceless modes, in particular. These sectors, however, are not involved within the 
minisuperspace sector of a homogeneous FRW metric, and it is suggestive to use a simple method of 
this finite ghost-avoidance renormalization which ultimately eradicates ghosts in this sector. 
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in the integrand of this expression will play a crucial role in the dynamical suppression 
of the vacuum Hartle-Hawking type instantons. 

5. Effective action on a static Einstein instanton 

Now we need r^fg] on a static background (TT7)|) with a periodically identified time 
1] having the period r/ . This is a typical thermodynamical calculation which, for 
completeness, we fully present here for a conformal scalar field. By decomposing it 
in spherical harmonics (enumerated by a collective index n) on a unit 3-sphere one gets 
the classical action as a sum of (Euclidean) oscillators with energies io n 

s[m = \^[<in({^)\^l). (25) 

The value of this action on solutions of classical equations 

^-u*<f> n = Q, n = l,2,..., (26) 

0n(O) = Vni 0n(%) = <^n, (27) 

gives the Hamilton- Jacobi function 

a/ I / n\ On + <Pn) COsh^p) - 2 <p n <p' n 

Sfa rj \(p,0) = 2^ "» 2sinh(u^ ) " (2§) 

Via the Pauli-Van Vleck formula it yields, up to a numerical factor, the evolution 
operator in the Euclidean time rjo - the density matrix 

D[0]exp(-S[S,0]) 

<t>n{m) = L Pn, <t>n(0) = (p' n 

Its trace generates the needed effective action 

e -rVi oop [s] = J Y[ d(fn J D[(j)] eX p(-S'[^0]) 



4>n(rio) = 4>n(0) = ipn 



= const Y\ (sinh j (3Q) 

n 

so that up to an additive constant it equals the sum of contributions of the vacuum 
energy E and free energy F(r} ) 

rVioopts] =J2 Y + ln(l-e-^ )] =m 2 P E oVo + F(r ]o ), (31) 

n 

oo 3 

^o^rE? (32) 

n n=l 

oo 

F(?7o) = J2 ln i 1 ~ e ~" nm ) = ^n 2 \n{l- e^™ 70 ). (33) 

n n=l 
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Here we go over to the explicit summation over the principal quantum number n = 
1, 2, 3, in terms of which u n = n and n 2 is the degeneracy of the energy level. Also 
we note that our definition of the free energy F(j]o) differs from the conventional one by 
a factor of r) . 

Quartic divergence of the vacuum energy (|3^|) in Fi_i oop [g] constitutes the 
ultraviolet divergences of the full action r^ioopfg]. Under a covariant regularization 
the power and quartic divergences among them are absorbed by the renormalization of 
the cosmological and Einstein terms, whereas the subtraction of logarithmic divergences 
yields as a remnant the contribution of a conformal anomaly considered above. For 
conformal fields, which we consider, the logarithmic divergences are actually zero, 
because they are given by the sum of integrated Weyl-squared and Euler number terms 
(7 and (3 terms of the conformal anomaly |T7jl ). For a conformally flat metric with 
the torus topology they are both vanishing. Therefore, our regularized one-loop action 
actually does not have a typical renormalization ambiguity quadratic in the curvature 
- the term of the same structure as EqT]q i n (j31|) . J d i xg l l 2 R 2 ~ J dra 3 /a 4: . Thus, the 
vacuum energy in an Einstein static spacetime (}3*2*j) should be uniquely calculable. This 
was independently confirmed by different methods [SI], and it equals 

1 f 4 

m2f,£o= 960 X 17 (34) 
I 88 

respectively for scalar, spinor and vector fields, It should be emphasized that (fHlj) 
is the renormalized vacuum energy obtained from (|3*2*|) by subtracting the quartic 
and quadratic divergences. The result of this renormalization is a positive Casimir 
energy (which is positive also for the naively negative vacuum energy of a spinor field 

-EnK/2), EH). 

Finally, we have to take into account the effect of the finite ghost-avoidance 
renormalization (|2^j) which should be applied also to r^oopfg]. Since J d^xg 1 / 2 ^ 2 = 
72 7r 2 77o, this renormalization leads to the expression similar to (J31j) . but with the 
modified vacuum energy which we denote by Co 

T R [g]=m 2 P C oVo + F(r ]o ), (35) 
3 

m 2 P Co = m 2 P E H a. (36) 

16 

One can directly check by using (J 19)) and ()34|) that for all conformal fields of low 
spins the modified energy reduces to the one half of the coefficient B in the conformal 

% This is the second major point of departure from [7||H]. Covariant regularization of the vacuum energy 
never generates terms of the form (M^/H 4 ) f dr/a in which the structure f dr/a is accompanied by 
powers of ultraviolet and infrared H 2 cutoffs advocated in IH] - This structure arises for generic 
non-conformal fields from logarithmically divergent curvature-squared terms only in the combination 
ln(M^ / H 2 ) f dr/a, because powers of the ultraviolet cutoff can only renormalize effective A (f dra 3 
structure) and G (f dr a structure) and do not accompany radiation-type terms. But logarithms are 
too weak in the infrared to compete with the tree- level exponential of l|T|l. 
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part of the total effective action (|24j) + 

m 2 P C = -m 2 P B. (37) 

This relation is confirmed by ^-functional and point-separation regularizations, but not 
the dimensional one. The latter results in the alternative value of a for a vector field, 
a = —2/15 instead of 1/5 in (|19j). and also gives different results for the renormalized 
vacuum energy. This lack of universality casts certain doubt on the dimensional 
regularization, and in what follows we will accept the values consistently obtained 
within other regularization schemes. They maintain the relation (|37|) which will be 
very important below. 



6. Effective Friedmann and bootstrap equations 

Now we assemble together the classical part of the action, contributions of the conformal 
anomaly (|23j) and of the static instanton (|35|). By rewriting the conformal time as a 
parametrization invariant integral in terms of the lapse N and the scale factor a, 

f T+ air N(t) , N 

Vo = 2 / —f-Y-, 38 

Jr. «( r ) 

we finally have the full effective action 

r[o(r), N{t) } = 2mp £ dr (-^ -Na + NH 2 a 3 ^j 

ir> 2 r+ fa 2 1 d 4 
+ 2Bm 2 P I dr[ — 



\Na 6 N 3 a 
7 drN\ „ 2 r+drX 



+ F 2 / + Bmp / . (39) 

V Jr_ a J Jr. a 

It is essentially nonlocal due to the nonlinear dependence of the free energy term on the 
conformal time ([38)1 . which incorporates the "breathing" of the instanton in r-direction 
caused by local variations in a(r) and N(r). 

The effective Friedmann equation obtained from this action by varying the lapse 
reads as 

5T n ( aa 2 



__ 2l4 (__ + „V 



n f a 2 1 d 4 \ 2 fdF(ri ) B 9 \ , in . 

+ 2Bm 2 P — + — + - r + —ml 440) 

p \ N 2 a 2 N 4 aJ a \ dr] 2 P J ^ ' 

One can easily check that the variation of a does not result in the new equation because 
of the one-dimensional covariance of the action amounting to the Noether identity 

d(r) 5T d 5T 

N{r) 5a{r) = d^5N{r)' ( ' 

+ This result implies that the vacuum energy of conformal fields in a static Einstein universe can be 
universally expressed in terms of the coefficients of the conformal anomaly m 2 P Eq = 3(2/3 — a)/16. 
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In the gauge N = 1 the effective equation reduces to the Friedmann equation (0) 
modified by the quantum B-term 

a 2 /la 4 d 2 \ _ 1 C 

^ + B \2^-^)-^- H (42) 

in which the radiation constant C is a functional of a(r), determined by the bootstrap 
equation 



m^ = ^| + ^M. (43) 
2 c/?7q 



Here 



dF{rjo) n 



3 



o 



y (44) 

n=l 

is the thermal energy of a hot gas of particles, which adds to their vacuum energy 
m 2 P B/2. Thus, the overall back reaction effect is mediated by the contribution of the 
radiation-type energy density term and the anomalous quantum I?-term. In view of Eqs. 
f!38|) and (}4*3*|) the constant C characterizing the amount of radiation nonlocally depends 
on the FRW background supported by the radiation itself, and this is the mechanism 
of the bootstrap we are going to analyze.* 

The on-shell action T on solutions of this bootstrap system can be cast into a 
convenient form by expressing the combination —a+H 2 a 3 in (|39|) in terms of other pieces 
of the effective Friedmann equation (}4~2*|) . Then, after the conversion of the integral over 
t into the integral over a between the turning points a± (if any) the on-shell action 
takes the form 

r„ = F( m ) - „ ^M + 4 m | ** ( B - a 2 - *f) . (45) 

The structure of the integral term here will be of crucial importance for the elimination 
of the infrared catastrophe and formation of the bounded cosmological landscape. 

7. The effect of conformal anomaly and bootstrap 

The modified Friedmann equation P2*|) can be solved for a 2 as 



,{a 2 - B) 2 2H 2 . 2 (a 2 -B) 
V B 2 ~B~ ^ ~ ^ ~ Q -) ~ B ■ (46) 

Therefore, it has the same two turning points a± as in the classical case (jHJ) providedjj 
a 2 _ > B, (47) 

* Note that, contrary to a conventional wisdom, the vacuum energy of quantum helds in a static 
universe l|84|l (and its renormalized version <|rtf>|) 'I contributes not to the cosmological constant, but to 
the total radiation density. This follows from the structure of the corresponding terms in the effective 
action discussed in Sect. 5, cf. footnoteEl and is also confirmed by the radiation-type equation of state 
derived for Eq in \M\ . 

t) Only the positive sign of the square root is admissible in 146JI because otherwise the Euclidean 
configuration exists only outside of the domain a_ < a < a+, which violates the periodicity of the 
instanton solution. 
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because in the opposite case the second term of (J46|) at a_ is positive and a(r) at the 
contraction phase without obstacle crosses the value a_. The scale factor continues 
decreasing and either becomes complex, at a = [B(2C — B)/(l — 2BH 2 )] 1 / 4 , or reaches 
zero. In the first case the instanton solution is not periodic and should be excluded, 
whereas the second case exactly corresponds to the formation of the pure-state instanton 
of the Hartle-Hawking type. Indeed, the full conformal time tends to infinity, i] — > oo, 
because of the divergence of the integral (JHHj) at a(r) — > a(r_) = 0, so that both 
F(Vo) ~ — exp(— rjo) and dF(r]o)/dr]o vanish. Therefore, from the bootstrap equation 
C = B/2, and the instanton smoothly closes at a = 0, because in view of (J4*B*)) 

d 2 -> 1 + y/1 - 2C/B = 1, a -> 0. (48) 

If it were a conical singularity, the instanton solution could have been ruled 
out for reasons similar to those discarding the Hawking- Turok instantons F° r 
smooth instantons this argument does not apply. Nevertheless, smooth Hartle-Hawking 
instantons arising outside of the domain (|4"7jl are ruled out by the infinitely large positive 
value of their action. For their solutions the action T reduces to the last integral term of 
(|4"5Jl with the lower limit a(r_) = 0. Due to the contribution of the conformal anomaly 
and in view of (J4*%|) its integrand is positive at a — > 0, and the integral diverges at the 
lower limit to +oo. This is the mechanism of how thermal fluctuations destroy the 
Hartle-Hawking pure state and make a mixed state of the Universe dynamically more 
preferable. 

Together with the upper bound on CH 2 Q, the condition (|4"7j) is equivalent to the 
following inequalities 

ACH 2 < 1, C > B - B 2 H 2 , BH 2 < -, (49) 

which automatically imply that C > B/2. Thus, the admissible domain for instantons 
reduces to the curvilinear wedge of the (if 2 , C)-plane below the upper hyperbolic 
boundary and above the lower straight line boundary on Fig0] This domain lies to 
the left of the critical point C = B/2, H 2 = 1/2B, at which these two boundaries touch 
and form a cusp. 

Inside this domain our bootstrap eliminates the infrared catastrophe of H 2 — > 0. 
Indeed, the conformal time ()38|) can be rewritten as an integral over a between the 
turning points a±. In terms of the rescaled integration variable x = Ha it reads as 

r x + d x 

Jx ~ xU&^P- + l(x 2 + - x 2 ){x 2 - xl) - ^ ' 



x ± = ^=^l±^l^q, (51) 

and, in fact, turns out to be the function of two dimensionless parameters q and p, 

q = ACH 2 , p = BH 2 . (52) 




Figure 4. The instanton domain in the (H 2 , C)-plane is located between bold 
segments of the upper hyperbolic boundary and lower straight line boundary. 
The first one-parameter family of instantons is labeled by k = 1. Families of 
garlands are qualitatively shown for k = 2,3,4. (1/2B, B/2) is the critical 
point of accumulation of the infinite sequence of garland families. 



For H 2 — > (q —>■ 0) it diverges to infinity at least logarithmically at the lower 
integration limitft x_ — > 0, 770 — > 00. Then it follows from the bootstrap equation 
(|33jl that C ^ B/2, which is impossible because in view of (JH1J) C > 5 at # 2 = 0. 
Thus the one-parameter family of instanton solutions in the admissible wedge of the 
(H 2 , C)-plane never hits the C-axes and can only interpolate between the points on the 
boundaries of the domain (J49j) with positive values of H 2 . 

For a conformal scalar field the numerical analysis gives such a family starting from 
the point at the lower line boundary with the parameters 

H 2 « 2.997m 2 ,, C « 0.004 m~ 2 , r « -0.1559, (53) 

and terminating at the point on the upper hyperbolic boundary 

H 2 « 12.968 m|, C w 0.0193 m p 2 , r ^ -0.0883. (54) 

The action r changes in this family non-monotonically - it has the maximum 
To — —0.063 approximately at 7]q = 5.03, so that maximally probable are the end- 

ff In terms of the integral over a this is the divergence at a+ ~ 1/H — * 00, because of the singularity 
of the rescaling x — Ha at H — * 0. Thus, as it should be, this is the infrared limit of infinitely large 
universes. 
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point configurations, though the variation in their probability weights is qualitatively 
negligible. 

The instanton (|54jl describes the creation of a static Einstein Universe of constant 
size, a = a± = l/(y/2H), filled by a hot gas of a conformally invariant field in an 
equilibrium state with the temperature 
1 H 

T = = — - 55 

ar] vrVl - 2BH 2 

which is analytically calculable as a function of H, because in the limiting case of 
ACH 2 = q — > 1, a_ — > a+, the conformal time integral (|50|) yields 

/ 2 (1-2H 2 B) / =V2(1-2^)tt. (56) 

Remarkably, for a static Universe the Euclidean action To reduces to the first two 
terms of (|45|) which form the Legendre transform of F(j]q). The entropy of the full system 
which includes coupled radiation and gravity, is given by the Legendre transformation 
of r with respect to the inverse temperature j3 = r = l/ar] . Therefore, it equals 
S = —F(r)o). Apparently, this result can be reconciled with the known relation between 
the free energy F (it differs by the factor of /3 from the conventional definition), energy 
E and entropy, F = (3E — S, by noting that the total energy of a spatially closed system 
is zero, E — 0. 

The above results for a scalar field with m 2 p -B sca i ar = 1/240 can also be obtained for 
other spins. We present them for a vector field with a much bigger value of the constant 
m 2 P -B V ector = 31/120. The lower and upper bounds on the instanton family then read 
respectively as 

H 2 w 1.063 m 2 P , C w 0.187m P 2 , T ^ -0.2797, (57) 
H 2 w 1.229 m 2 P , C w 0.203 m P 2 , T w -0.2334. (58) 

Together with (|53|) they give the lower bounds on the instanton range for scalars and 
vectors ©-© cited in Introduction. Their comparison shows the tendency of decreasing 
Hubble scale and growing — ro with the growth of spin (or the parameter B). 

The corresponding behavior of our bootstrap at B ^> 1 can be studied analytically. 
For this purpose it is enough to stay at the hyperbolic boundary, H 2 = 1/4(7, where 
the value of the conformal time is given by ()56|) . In view of the bootstrap equation 
it equals 



. F'( m ) _ dF(r h ) 

ri0 = 2 ^ m ^B + 2F'( Va y F( » ) = ^T' (59) 

For growing B the solution of this equation tends to zero, 770 — > 0, so that one can 
approximate the high-temperature statistical sums and (j4^|) by F(r] ) ~ — 7t 4 /45?7q 
and F'(j]o) ~ 7t 4 /15?7q. Therefore, from (jBTIj) one has the asymptotic behavior of the 
conformal time r] ~ 7r(4/15 m 2 P B) 1 ^ really decreasing with the growth of B. For static 
instantons at the upper boundary the integral term of ()45|) identically vanishes, and the 
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on-shell action reduces to the contribution of the thermal terms 

T ~ = m P VB, m 2 P B » 1. (60) 

3v 15 

It is important that To stays negative (even though no infrared catastrophe exists any 
more) while its absolute value grows with B. This property implies that the growing 
spin of quantum matter (growing B = 3(3 /imp, cf. ()19|)) makes the probability of the 
underlying instanton higher. Also it will be important for suppressing the contribution 
of Lorentzian (overbarrier) configurations in extended formulations of quantum gravity, 
briefly discussed in Conclusions. 

Qualitatively, this behavior stays true also for instantons below the upper boundary. 
Indeed, with a growing B the lower straight line boundary of ()49j) becomes steeper, the 
critical point (1/2B, B/2) slides along the hyperbolic boundary in the upward direction. 
The whole family of bootstrap solutions drops deeper into the curvilinear wedge of Fig0] 
approaching closer to this point, (C — B/2) /(B/2) ~ (30-B) -1 / 3 — > 0, in the vicinity of 
which members of this family do not differ much from one another (cf. Eqs. (l57jl - ()58jl ). 

The limit of large B is very important. As will be discussed in Conclusions, this 
limit arises when ascending the hierarchy of particle phenomenology with a growing spin 
and number of matter multiplets. This limit corresponds to strong quantum corrections 
mediated by a large contribution of the conformal anomaly. These corrections lead to 
a rather unexpected result. The higher the value of B, the stronger it truncates the 
instanton domain (|4"9*j) from above, H 2 < -f^ax = 1/25. This property will be discussed 
in Conclusions as a possible mechanism of reducing the landscape of string vacua. 

8. Instanton garlands 

The instanton family of Sect. 7 does not exhaust the entire cosmological landscape. 
Together with simple instantons of the above type one can consider multiple folded 
instantons in which the scale factor reaches its maximal and minimal values many times, 
similarly to the instantons considered in [33] ■ They can be obtained by glueing together 
into a torus (at surfaces of the maximal size a + ) k copies of a simple instanton, see 
FigE3 We will call them garlands. Their sequence will generate a new quantum scale 
in gravity theory — the upper bound of the cosmological landscape (J2J). 




Figure 5. Segment of the garland consisting of three folds of a simple instanton 
glued at surfaces of a maximal scale factor. 
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The full lapse of the Euclidean conformal time for a fc-length garland equals 

(k) , f T+ dr , f a + da , , 

4 k) = 2k — = 2k — , 61 

Jr. a Ja. aa 

and the whole bootstrap formalism persists with the replacement of ()38|) by rj^ and 
multiplying the last (integral) term of ()45|) by k. The numerical analysis for k = 2 
shows the existence of the instanton family similar to the case of k = 1. It interpolates 
between the point on the lower boundary of (C, if 2 )-plane 

Hf 2) « 45.89 m%, C m « 0.0034 nip 2 , T {2) « -0.0113, (62) 

and the point on the upper (hyperbolic) boundary 

# ( 2 2) « 61.12 m 2 ,, C (2 ) « 0.0041 mp 2 , r{, 2) « -0.0145. (63) 

Moreover, it turns out that such families exist for all 1 < k < oo, and their 
infinite sequence tends to the critical point C = B/2, H 2 = 1/2B, where the lower and 
upper boundaries merge. For each k the instanton family forms a line joining the upper 
and lower boundaries, the length of this line getting shorter and shorter with k — > oo 
and its location closer and closer approaching the critical cusp point. As in the case 
of k = 1 the upper point of each family gives rise to a hot static Universe filled by 
radiation in the equilibrium state with the temperature T^) differing from (|55|) by an 
extra 1/k- factor. 

The existence of this sequence of instanton families follows from the behavior of 
the single fold of the conformal time ()56|) at 2BH 2 — > 1. It tends to zero in this limit, 
so that when multiplied by k it admits the solution of the bootstrap equation for any 
k — > oo with the total time r]^ slowly growing to infinity. In this limit the sequence 
of instanton families can be described analytically. For simplicity consider the upper 
hyperbolic boundary where e = C — B/2 tends to zero as we approach a critical point 
and H 2 ~ (1 — 2e)/2B. The full conformal time (that is ()56|) times k) then behaves as 
7^ ~ 2nk^l ej B. In the assumption that 77^ — * 00 for k — * 00 the bootstrap equation 
reads as 

uipE ~ exp(— 1]^), (64) 

because in this limit F'(t]q) ~ exp(— n ). It has the solution, e ~ Sin 2 /c 2 /47r 2 fc 2 , which 
indeed confirms that the conformal time is slowly growing to infinity, r]^ ~ \nk 2 — > 00. 
A similar analysis can be performed for the lower boundary of the instantin domain. 
Thus, it turns out that within the l//c 2 -accuracy both upper and lower points of each 
k-th family coincide and read as 

2 1 / ln 2 k 2 \ 
^M 1_ 2W> (65) 
B( ln 2 fc 2 \ , , 

rf 1 -^^. (6T) 
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Thus, the lengths of instanton families (both in H 2 and C directions) decrease as the 
second order of the l//c 2 -expansion, 1/k 4 , so that on Fig0]they fit in ever narrowing 
wedge near the critical cusp of C = B/2 = 1/AH 2 . With a growing k all garlands 
become more and more static and cool down to zero temperature 

T «" = 7W^ - °- (68) 



It is remarkable, that in contrast to the tree-level instantons of [o3] the garland 
action is not additive in k, so that with k — > oo it goes to zero rather than to — oo. This 
happens because a naively additive part of the action - the integral term of ()45|) even 
when multiplied by k tends to zero with a growing number of folds. This in its turn 
follows from the growingly static nature of garlands - in terms of the deviation from the 
critical point e we have a + — a_ ~ e, a ~ yfe, B — a 2 ~ e, and the integral term behaves 
like k (lnfc 2 ) 5,/2 //c 4 . Because of the decreasing temperature ()68|) the thermal part 

of flj3)) also decreases, but slower than the integral term, and yields (JB7J). 

Thus, infinitely long garlands do not dominate the instanton distribution, but their 
existence is very important, because they generate a new quantum scale — the upper 
bound of the instanton range (J2J). Their sequence converges to the cold, T^) = 0, static 
instanton with the vanishing action = 0, which realizes this scale as a maximal 
possible value of the Hubble constant in the instanton landscape 

#Lx = A?*) = 1/25. (69) 

For a scalar and vector fields it gives respectively the values 120 m 2 P and 60/31 m 2 P cited 
in Introduction, ©-Q- 



9. Conclusions 



The picture of the cosmological landscape that we get from the back reaction of quantum 
matter is very attractive. It consists of the ensemble of universes in the limited range 
of the cosmological constant, filled by a hot quasi-equilibrium gas of conformal field 
particles. After nucleation from instantons these universes expand in the Lorentzian 
regime. Eventually the radiation dilutes to a negligible density when the cosmological 
constant A = 3H 2 starts dominating and supports the inflationary stage, after which 
we get a usual large-scale structure formation scenario. This mechanism works also for 
static Einstein universes (originating from the upper boundary of the instanton domain), 
because they also start expanding as a result of their gravitational instability. 

It is interesting that the Lorentzian version of the effective Friedmann equation 
(142)1 admits at late stages of evolution, a ^> 1/H , two branches of the quasi-exponential 
expansion with the Hubble parameters 

H 2 ± = — (1± Vl -2BH 2 ). (70) 
B 

The original Hubble constant H should have a composite (inflaton-induced) nature 
in order to decay at the exit from inflation, so that if 2 ~ H 2 — > 0, whereas the 



Cosmological Landscape From Nothing 



22 



second branch enters the acceleration stage with — > 2/B. Unfortunately, for 
hierarchy reasons this mechanism cannot play the role of dark energy, because the 
needed value of B should be 10 120 nip 2 . It is tempting to speculate, though, that at the 
fundamental level the parameter B is a moduli variable which is sufficiently small in 
the early Universe, but becomes exponentially big at late stages due to some dynamical 
compactification mechanism of the braneworld type [31] (see also [35] for branes in 
Euclidean quantum gravity). This would imply that the new quantum gravitational 
scale, induced by garlands, is responsible for the dark energy phenomenon. However, 
the detailed consideration of this idea goes beyond the scope of this paper. 

Our cosmological landscape belongs to the bounded range (J2J. Its infrared cutoff 
is provided by the radiation back reaction which makes a mixed state of the early 
Universe dynamically more preferable. Back reaction dynamically suppresses the pure 
state Hartle-Hawking instanton by its infinitely large positive action generated by the 
conformal anomaly. In essence, this implies a solution to the long-standing problem 
of boundedness of the on-shell cosmological action in Euclidean quantum gravity and 
eliminates the infrared catastrophe of a small cosmological constant (This should not be 
confused with the off-shell indefiniteness of the Einstein action caused by the metric 
conformal mode). This effect survives even in the classical limit of the bootstrap 
equations ([42 )1 - (|43 |) . B = 0, for which the lower bound of (J2J can be obtained in an 
analytical form. In this limit the conformal time (|5()jl for q = 4CH 2 — ► behaves 
as 7/0 — — ln(g/64) — 3glng/4 — ► oo. When substituted into the bootstrap equation 
(@3J with F'(t/ ) ~ exp(-77 ), it gives H 2 ~ 16ml [1 + 48m P Cln(l/m P C)} -> 16m 2 P , 
C — > 0. Therefore, the infrared catastrophe is ruled out already at the classical level by 
the infrared cutoff which for a single conformal scalar field equals A m j n = 48 m 2 P . 

In contrast, the high-energy cutoff is the quantum effect of vacuum energy 
and conformal anomaly, which generates a new scale in gravity theory. In view of the 
relation ([21]) this scale is determined by the inverse of the coefficient of the topological 
Gauss-Bonnet term in the conformal anomaly (|17|). 

2 

A max = -m P , m 2 P = 3ir/AG, (71) 

and it tends to infinity in the classical limit (3 — > 0. In the cosmological landscape this 
scale is realized as a limiting point of the sequence of garland-type instantons. 

Of course, we considered only the contribution of conformally invariant fields of a 
few low spins. Other fields, gravitons in particular, can change the picture, but it is 
suggestive that the situation will qualitatively stay the same, because on static or quasi- 
static backgrounds their behavior does not differ much from their conformally invariant 
analogues. What is very critical for our results is the value of the vacuum energy E 
of quantum fields, which after the ghost-avoidance renormalization yields a particular 
contribution Cq = B/2 to the total constant C in the bootstrap equation (|43|). Quite 
interestingly, this value guarantees that the bootstrap solutions have a lower bound for 
C exactly coinciding with the critical point C = B/2 = 1/AH 2 at which the infinite 
sequence of garlands is accumulating. This peculiar fact apparently implies a sound link 
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between the renormalization theory on curved spacetime background, ghost avoidance 
criterion and quantum gravitational tunneling. 

Conformally non-invariant fields are likely to destroy this relation, and one can 
imagine three possibilities. If overall Co is less than overall B/2 the picture qualitatively 
stays the same. In particular, an infinite sequence of garlands survives, though it tends 
the critical point at a finite value of the conformal time rjo and, therefore, with a finite 
temperature. If B > C > B/2, this infinite sequence becomes truncated at some 
maximal number of instanton folds. Finally, if Co > B the infrared catastrophe of 
A — > occurs again - the one-parameter family of instantons necessarily hits the C- 
axes of the (C, iJ 2 )-plane at Co- This point exactly solves the bootstrap equation and 
infinitely dominates the rest of the partition function (at H 2 — > the upper integration 
limit in ()45)1 a + — > oo, and the dominant second term of the integrand at this limit 
makes To — > — oo). Which of these three possibilities is realized is an open question 
deserving further study. 

Another open question concerns the normalizability of our partition on the infinite 
set of instantons. One might think that it is not normalizable because of infinite 
summation over garland folds k. However, at least naively the total continuous measure 
of the instanton range (J2J) is finite because (Amix — A ffi n ) ~ 2~2k(^-/^ 4 ) < 00 • m 
order to have a definitive conclusion, though, one must take into account a zero-mode 
contribution to the actual measure and also estimate preexponential factors. The latter 
reduce to the quantum mechanical functional determinants of the Hessian of the effective 
action r[a(r), N(r) ], which is a nonlocal operator rather than a differential one. Their 
calculation is doable and will be reported elsewhere |3l)] . 

The boundaries of the instanton range (J2J) depend on particle phenomenology. For 
a single conformal or vector field they have Planckian values ©-©• Therefore, with the 
nonperturbative back reaction taken into account, we still have the issue of semiclassical 
validity of our model. This problem can be handled by adding new multiplets of quantum 
fields and raising their spin. There is a simple scaling behavior in the number of fields 
N. All terms in the bootstrap equation ()43)) rescale linearly in N, 

C^NC, B-^NB, F(rj ) -> NF(r]o), (72) 

whereas tjq given by (|5Uj) depends only on combinations BH 2 and CH 2 . Therefore, 
the bootstrap remains consistent when in addition to ()72|) the value of H 2 scales down 
as H 2 — > H 2 /N. This justifies the semiclassical approximation for large N with the 
relevant 1/A^-expansion. 

Moreover, as was noted in the end of Sects. 7 - 8, the growing spin of a conformal 
particle (or the parameter B = 3/3/4 mp) also lowers the landscape scale and in view of 

makes its probability weight higher. This means that when ascending the hierarchy 
of spins, this scaling in N and B reduces the domain (J2J) to more and more narrow 
subplanckian range and suggests a long-sought selection mechanism for the landscape 
of string vacua. Modulo the details of a relevant 4D-compactification, this mechanism 
might work as follows. For (3 rapidly growing with N, [3 ~ N, and with the spin, cf. Eq. 
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(Unj), the upper scale (|7Tj) decreases towards the increasing phenomenology scale, and 
converges to the latter at the string scale m 2 , 1/B — ► m 2 , where a positive A might be 
generated by the KKLT or KKLMMT-type mechanism [37] . Our conjecture is that at 
this scale our landscape of hot instantons selects from the enormous landscape of string 
vacua a reasonably small subset compatible with the observed particle phenomenology 
and large-scale structure. The corresponding bootstrap is perturbatively consistent 
provided 1/B = m 2 <^ m 2 P . 

Our results hold within Euclidean quantum gravity incarnated in the path integral 
(110}) over Euclidean configurations. It gives rise to the Lorentzian spacetime by 
analytical continuation but, otherwise, automatically excludes Lorentzian saddle point 
configurations. However, one can imagine an extended formulation of quantum gravity 
generalizing the path integral (fT0|) to a wider integration domain or using an alternative 
formalism, like the Wheeler-DeWitt equation. Within such an extension Lorentzian 
solutions can, in principle, exist above the upper hyperbolic boundary of (J4U)) . in the 
domain 4CH 2 > 1 overbarrier from the viewpoint of the minisuperspace Wheeler- 
DeWitt equation. They will contribute to the cosmological density matrix, and this 
casts certain doubt on the universality of our results. Our conclusions, however, are 
likely to stay true. Indeed, according to ()60|) for a high value of B the effective action 
scales as r ~— mp\f~B. Therefore, because it is still negative our Euclidean landscape 
at the string scale mentioned above is weighted by exp(#m P ^/B) = exp(# mp/m s ) ^> 1. 
Therefore, it exponentially dominates over Lorentzian configurations, the amplitudes of 
which being 0(1) in view of their pure phase nature. Thus, our results look robust 
against possible generalizations of Euclidean quantum gravity. 

This is how cosmological landscape emerges from nothing and perhaps tames its 
string counterpart, provided some like it hot. 
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